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Abstract 

We study the quark deconfinement phase transition in hot /3-stable hadronic matter. Assuming a first order phase 
transition, we calculate the enthalpy per baryon of the hadron-quark phase transition. We calculate and compare 
the nucleation rate and the nucleation time due to thermal and quantum nucleation mechanisms. We compute 
the crossover temperature above which thermal nucleation dominates the finite temperature quantum nucleation 
mechanism. We next discuss the consequences for the physics of proto-neutron stars. We introduce the concept of 
limiting conversion temperature and critical mass M cr for proto-hadronic stars, and we show that proto-hadronic 
stars with a mass M < M cr could survive the early stages of their evolution without decaying to a quark star. 
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1. Introduction 

In the last few years there has been a growing 
interest in the study of the nucleation process 
of quark matter (QM) in the core of massive 
neutron stars. In particular, it has been shown 
1 2 3 4 5 6 7 8 that above a threshold value of 
the central pressure a pure hadronic compact 
star (HS) is metastable to the decay (conver- 
sion) to a quark star (QS) (i.e. to a hybrid neu- 
tron star or to a strange star [9110] . depending 
on the details of the equation of state (EOS) for 
quark matter used to model the phase transi- 
tion [11112113114] ). This stellar conversion pro- 
cess liberates a huge amount of energy (a few 
10 53 erg) [15] and it could be the energy source 
of some of the long Gamma Ray Bursts (GRBs). 

The research reported in Refs 1 2 3 4 5 6 7 8 
has focused on the quark deconfinement phase 
transition in cold (T = 0) and neutrino-free 
neutron stars. In this case the formation of the 
first drop of QM could take place solely via a 
quantum nucleation process. 

A neutron star at birth (proto-neutron star) 



is very hot (T = 10 - 30 MeV) with neutri- 
nos being still trapped in the stellar interior 
[16117] . Subsequent neutrino diffusion causes 
deleptonization and heats the stellar matter to 
an approximately uniform entropy per baryon 
S =1 — 2 (in units of the Boltzmann's constant 
&b). Depending on the stellar composition, 
during this stage neutrino escape can lead the 
more "massive" stellar configurations to the 
formation of a black hole |18|17] , However, if 
the mass of the star is sufficiently small, the star 
will remain stable and it will cool to tempera- 
tures well below 1 MeV within a cooling time 
tcooi ~ a few 10 2 s, as the neutrinos continue 
to carry energy away from the stellar material 
[16|17] , Thus in a proto-neutron star, the quark 
deconfinement phase transition will be likely 
triggered by a thermal nucleation process. In 
fact, for sufficiently high temperatures, thermal 
nucleation is a much more efficient process with 
respect to the quantum nucleation mechanism. 

Some of the earlier studies on quark matter 
nucleation (see e.g., [19I20I21I22I2"3"] ) have al- 
ready dealt with thermal nucleation in hot and 
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dense hadronic matter. In these studies, it was 
found that the prompt formation of a critical 
size drop of quark matter via thermal activation 
is possible above a temperature of about 2 — 3 
MeV. As a consequence, it was inferred that 
pure hadronic stars are converted to quark stars 
within the first seconds after their birth. How- 
ever, these works |19I20I21I22| reported an esti- 
mate of the thermal nucleation based on "typ- 
ical" values for the thermodynamic properties 
characterizing the central part of neutron stars. 

Our main objective in this and next related 
works, is to establish if a newborn hadronic star 
(proto-hadronic star) could survive the early 
stages of its evolution without " decaying" to a 
quark star. In the present Letter, we calculate 
the thermal nucleation rate of quark matter in 
hot (T ^ 0) and neutrino-free hadronic matter 
using a finite temperature EOS for hadronic 
and quark matter. In addition, we calculate 
the quantum nucleation rate at finite temper- 
ature, and compare the thermal and quantum 
nucleation time at different temperatures and 
pressures characterizing the central conditions 
of metastable proto-hadronic compact stars. 
We compute the crossover temperature above 
which thermal nucleation dominates above the 
finite temperature quantum nucleation mech- 
anism. Finally we briefly discuss some conse- 
quences for the physics of proto-neutron stars. 

2. Phase equilibrium 

For a first-order phase transition the con- 
ditions for phase equilibrium are given by the 
Gibbs' phase rule 

Tr = Tq = T , P H = P Q =P (1) 
Hh(T,P )=»q{T,P ) (2) 
where 

£h + Ph - s H T e Q + P Q - s Q T 
— , MQ — 

are the Gibbs' energies per baryon (average 
chemical potentials) for the hadron and quark 
phase respectively, e H (sq), Ph (Pq), s h (sq) 
and uh (no) denote respectively the total {i.e., 
including leptonic contributions) energy den- 
sity, total pressure, total entropy density, and 

1 We assume the quark deconfinement phase transition 
to be of the first order. This assumption is common in 
most of the studies of quark deconfinement in compact 
stars. 



baryon number density for the hadron (quark) 
phase. Above the "transition point" (Jo) the 
hadronic phase is metastable, and the stable 
quark phase will appear as a results of a nucle- 
ation process. 

Small localized fluctuations in the state vari- 
ables of the metastable hadronic phase will give 
rise to virtual drops of the stable quark phase. 
These fluctuation are characterized by a time 
scale Vq 1 ~ 10~ 23 s. This time scale is set by the 
strong interactions (which are responsible for 
the deconfinement phase transition), and it is 
many orders of magnitude shorter than the typ- 
ical time scale for the weak interactions. There- 
fore quark flavor must be conserved during the 
deconfinement transition. We will refer to this 
form of deconfincd matter, in which the flavor 
content is equal to that of the /3-stable hadronic 
system at the same pressure and temperature, 
as the Q*-phase. Soon afterward a critical size 
drop of quark matter is formed, the weak in- 
teractions will have enough time to act, chang- 
ing the quark flavor fraction of the deconfined 
droplet to lower its energy, and a droplet of j3- 
stable quark matter is formed (hereafter the Q- 
phase) . 

This first seed of quark matter will trigger 
the conversion [24125115] of the pure hadronic 
star to a hybrid star or to a strange star. Thus, 
pure hadronic stars with values of the cen- 
tral pressure larger than Pq are metastable to 
the decay (conversion) to hybrid stars or to 
strange stars [112131415] , The mean lifetime of 
the metastable stellar configuration is related 
to the time needed to nucleate the first drop of 
quark matter in the stellar center and depends 
dramatically on the value of the stellar central 
pressure [l|2|3|4|5j . 



3. Quantum and thermal nucleation 
rates 

The main effect of finite temperature on the 
quantum nucleation mechanism of quark mat- 
ter is to modify the energy barrier separating 
the quark phase from the metastable hadronic 
phase. This energy barrier, which represents the 
difference in the free energy of the system with 
and without a Q*-matter droplet, can be writ- 
ten as 

U(H, T) = ~7rn Q . (jjlq, - fi H )Tl 3 + AnaTZ 2 (4) 



2 



where 1Z is the radius of the droplet (supposed 
to be spherical) , and a is the surface tension for 
the surface separating the hadron from the Q*- 
phase. The energy barrier has a maximum at the 
critical radius 1Z C — 2a /[nQ*(nH — HQ*)]- No- 
tice that we have neglected the term associated 
with the curvature energy, and also the terms 
connected with the electrostatic energy, since 
they are known to introduce small corrections 
[29)2] . The value of the surface tension a for 
the interface separating the quark and hadron 
phase is poorly known, and typically values used 
in the literature range within 10 — 50 MeV fm -2 
[26128129] . We assume a to be temperature in- 
dependent and we take a = 30 MeV fm -2 . 

The quantum nucleation time r q can be 
straightforwardly evaluated within a semi- 
classical approach [27 28 29 . First one com- 
putes, in the WKB approximation, the ground 
state energy Eq and the oscillation frequency 
of the drop in the potential well U (1Z, T) . 
Then, the probability of tunneling is given by 



p = exp 



A(Eo) 



(5) 



where A(E) is the action under the potential 
barrier, which in a relativistic framework reads 



A(E) = - J V[2m(K)c 2 +E- U(JZ)][U{TZ) - 

TZ- 

being 1Z± the classical turning points and m(lZ) 
the droplet effective mass. The quantum nucle- 
ation time is then equal to 



(voPoNc)- 1 



(7) 



with N c ~ 10 48 being the number of nucle- 
ation centers expected in the innermost part 
(r < R nuc ~ 100 m) of the hadronic star, where 
the pressure and temperature can be considered 
constant and equal to their central values. 

The thermal nucleation rate can be written 
[321 as 



/=— Sl eM-U(Tl c ,T)/T) 



(8) 



where k is the so-called dynamical prefactor, 
which is related to the growth rate of the drop 
radius TZ near the critical radius (1Z C ), f^o is the 
so-called statistical prefactor, which measures 
the phase-space volume of the saddle-point re- 
gion around TZ C , and U(1Z C ,T) is the activa- 
tion energy, i.e. the change in the free energy 



of the system required to activate the forma- 
tion of a critical size droplet. The Langer the- 
ory [301311 32 33 of homogeneous nucleation has 
been extended in Refs. [34135] to the case of first 
order phase transitions occurring in relativistic 
systems, as in the case of the quark deconfinc- 
ment transition. The statistical prefactor, can 
be written [34] as 



fin 



^ a \^/2 ^ 



3V3VTV \ti Q J 



(9) 



where £q is the quark correlation length, which 
gives a measure of the thickness of the inter- 
face layer between the two phases (the droplet 
"surface thickness"). In the present calculation 
we take £q = 0.7 fm according to the estimate 
given in Refs. [34122] , 

For the dynamical prefactor we have used a 
general expression which has been derived by 
Venugopalan and Vischer [55] (see also Refs. 

ran 



2(7 



K 3 c {Aw) 2 



(10) 



where Aw = wq*—wh is the difference between 
the enthalpy density of the two phases, A the 
thermal conductivity, r\ and £ are the shear and 
bulk viscosities respectively of hadronic mat- 

ter. Notice that the nucleation prefactor used in 

E] (&he present work differs significantly from the 
one used in previous works [19120121] ) where, 
based on dimensional grounds, the prefactor 
was taken to be equal to T 4 . 

There are not many calculations of the trans- 
port properties of dense hadronic matter. With 
a few exceptions (see e.g. [37138] ). most of them 
are relative to nuclear or pure neutron matter 
[39140141142155] . These quantities have been cal- 
culated by Danielewicz [30] in the case of nu- 
clear matter. According to the results of Ref. 
[40] , the dominant contribution to the prefactor 
k comes from the shear viscosity r). Therefore, 
we take A and C equal to zero, and we use for 
the shear viscosity the following relation [40] : 



V 



7.6 x 10 26 



(T/MeV) 2 V n 



MeV 
fms 



(11) 



with no = 0.16 fm -3 being the saturation den- 
sity of normal nuclear matter. 

The thermal nucleation time r t h, rela- 
tive to the innermost stellar region (V nuc = 
(47r/3)i? 3 iuc ) where almost constant pressure 
and temperature occur, can thus be written as 
rth = (Vnuciy 1 - 
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Fig. 1. (Color online.) Gibbs energy per baryon as a 
function of pressure for the hadronic and quark phases 
at different temperatures. Lines with the larger slope 
refer to the hadronic phase. Full dots indicate the tran- 
sition pressure Pq for each temperature. 




Fig. 2. Phase equilibrium curve for the hadron to quark 
matter phase transition. 

5. Results and discussion 



4. Equation of state 

Over the last decade, it has been realized that 
strong interacting matter at high density and 
low temperature may possess a large assortment 
of phases. Different possible patterns for color 
superconductivity have been conjectured (see 
e.g. [44 45J and references therein quoted). Very 
recently, a new phase of QCD, named quarky- 
onic phase, has been predicted |46|47j . This hy- 
pothetical matter phase is characterized by chi- 
ral symmetry and confinement [46 47 48 . 

In the present work, we have adopted a 
more traditional view, assuming a single first 
order phase transition between the confined 
(hadronic) and deconfined phase of dense mat- 
ter, and we have used rather common models 
for describing them. For the hadronic phase we 
have used models which are based on a relativis- 
tic Lagrangian of hadrons interacting via the 
exchange of cr, p and ui mesons. We have used 
one of the parameters sets given in Ref. [49] : 
hereafter we refer to this model as the GM1 
equation of state. For the quark phase we have 
adopted a phenomenological EOS [H] which 
is based on the MIT bag model for hadrons. 
In this work, we have used the following set of 
parameters: m u = m ( i — 0, m s = 150 MeV for 
the masses of the up, down and strange quark 
respectively, B = 85 MeV/fm 3 for the bag con- 
stant, and a s = for the QCD structure con- 
stant. The two models for the EOS have been 
generalized to the case of finite temperature. 



In Fig. 1 we plot the Gibbs' energies per 
baryon for the hadron-phase and the Q*-phase 
at different temperatures, T = 0, 10, 20, 30 
MeV; lines with the larger slope refer to the 
hadron-phase. As we see, the transition pres- 
sure Po (indicated by a full dot) decreases 
when the hadronic matter temperature is in- 
creased. The phase equilibrium curve Po(T) 
for the hadron-quark phase transition (within 
the present schematic model for the EOS) is 
shown in Fig. 2. As it is well known, for a first- 
order phase transition the derivative dPo/dT 
is related to the specific (i.e. per baryon) la- 
tent heat Q of the phase transition by the 
Clapeyron-Clausius equation 

dP = n H n Q , Q 
dT iiq* — nn T 

Q = W Q *-W H = T(S Q * ~ S H ) (13) 

where Wh (Wq*) and Sh (Sq*) denote respec- 
tively the enthalpy per baryon and entropy per 
baryon for the hadron (quark) phase. The spe- 
cific latent heat Q and the phase numbers den- 
sities nn and uq* at phase equilibrium are re- 
ported in Table 1. As expected for a first order 
phase transition one has a discontinuity jump 
in the phase number densities: in our particular 
case nQ*(T,P ) > uh(T,P ). This result, to- 
gether with the positive value of Q (i.e. the de- 
confinement phase transition absorbs heat) tell 
us (see Eq. (fl"2j) ). that the transition tempera- 
ture decreases with pressure (as in the melting 
of ice). 
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0.358 


38.58 


15 


5.71 


0.439 


0.348 


36.55 


20 
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0.334 


33.77 
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17.17 


0.414 


0.316 


30.36 


30 


25.44 


0.398 


0.294 


26.53 


Table 1 



The specific latent heat Q and the phase numbers den- 
sities rifi and tiq* at phase equilibrium. 
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Fig. 3. Energy barrier for a virtual drop of the Q*-phasc 
as a function of the droplet radius and for different 
temperatures for a pressure P = 57 MeV/fm 3 . 

In Fig. 3, we represent the energy barrier for 
a virtual drop of the Q*-phase as a function 
of the droplet radius and for different tempera- 
tures at a fixed pressure P — 57 MeV/fm 3 . As 
expected, from the results plotted in Fig. 1, the 
energy barrier U(TZ, T) and the droplet critical 
radius 1Z C decrease as the matter temperature 
is increased. This effect favors the Q*-phase for- 
mation, and in particular increases (decreases) 
the quantum nucleation rate (nucleation time 
T q ) with respect to the corresponding quantities 
calculated at T = 0. 

In Fig. 4 we plot the quantum and thermal 
nucleation times of the Q*-phase in /3-stable 
hadronic matter as a function of temperature 
and at a fixed pressure P = 57 MeV/fm 3 . As 
expected, we find a crossover temperature T co 
above which thermal nucleation is dominant 
with respect to the quantum nucleation mech- 
anism. For the case reported in Fig. 4, we have 
T co = 7.05 MeV and the corresponding nucle- 
ation time is log 10 (r/s) = 54.4. The crossover 




T [MeV] 

Fig. 4. (Color online.) Thermal (rth) and quantum (t 9 ) 
nucleation time of quark matter (Q*-phase) in /3-stable 
hadronic matter as a function of temperature at fixed 
pressure P = 57 MeV/fm 3 . The crossover tempera- 
ture is T co = 7.05 MeV. The limiting conversion tem- 
perature for the proto-hadronic star is, in this case, 
© = 10.3 MeV, obtained from the intersection of the 
thermal nucleation time curve (continuous line) and the 
dot-dashed line representing log 10 (r/s) = 3. 
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53.98 


5.0 


233.6 


55.48 


6.0 


121.3 


56.94 


7.0 


56.6 


58.42 


8.0 


16.0 


58.85 


8.3 


3.0 


Table 2 



Crossover temperature T co (in MeV), for different fixed 
values of the pressure P (in Mev/fm 3 ) of hadronic mat- 
ter. The third column reports the logarithm of the nu- 
cleation time (in seconds) calculated at the crossover 
temperature. The value 8.3 MeV defines the value of 
the limiting conversion temperature for a star with 
a central pressure P = 58.85 MeV fm — 3 . 

temperature, for different values of the pressure 
of /3-stable hadronic matter, is reported in Tab. 
2 (second column) together with the nucleation 
time calculated at T — T co (third column). 

Having in mind the physical conditions in 
the interior of a proto-hadronic star [16ll7j (see 
Sect. I of the present paper), to establish if this 
star will survive the early stages of its evolu- 
tion without "decaying" to a quark star, one 
has to compare the quark matter nucleation 
time t = mm.{T q ,Tth) with the cooling time 
tcooi ~ a few 10 2 s. If t >> tcoal then quark 
matter nucleation will not likely occur in the 
newly formed star, and this star will evolve to a 
cold dclcptonizcd configuration. We thus intro- 
duce the concept of limiting conversion temper- 
ature for the proto-hadronic star and define 
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Fig. 5. (Color online.) The limiting conversion temper- 
ature for a newborn hadronic star as a function of the 
central stellar pressure. Newborn hadronic stars with a 
central temperature and pressure located on the right 
side of the curve @(P) will nucleate a Q*-matter drop 
during the early stages of their evolution, and will fi- 
nally evolve to cold and deleptonized quark stars, or will 
collapse to black holes. The lines labeled Tg represent 
the stellar matter temperature as a function of pressure 
at fixed entropies per baryon S/kg = 1 (dashed line) 
and 2 (solid line). 

it as the value of the stellar central temperature 
T c for which the Q*-matter nucleation time is 
equal to 10 3 s. The limiting conversion temper- 
ature O will clearly depend on the value of the 
stellar central pressure (and thus on the value 
of the stellar mass). 

The limiting conversion temperature O is 
plotted in Fig. 5 as a function of the stellar 
central pressure. A proto-hadronic star with a 
central temperature T c > O will likely nucleate 
a Q*-matter drop during the early stages of its 
evolution, and will finally evolve to a cold and 
deleptonized quark star, or will collapse to a 
black hole (depending on the particular model 
adopted for the matter EOS). 

For an isoentropic stellar core |16|17j . the 
central temperature of the proto-hadronic star 
is given, for the present EOS model, by the 
lines labeled by T$ in Fig. 5, relative to the case 
S = 1 ks (dashed curve) and S = 2 ks (con- 
tinuous curve). The intersection point (Ps, 0s) 
between the two curves @(P) and T$(P) thus 
gives the central pressure and temperature of 
the configuration that we denote as the critical 
mass configuration of the proto-hadronic stellar 
sequence. The value of the gravitational critical 
mass M cr = M(Pg,Qs) and baryonic critical 
mass Mb,ct are reported in Tab. 3, for three 
different choices of the entropy per baryon, 



S/k B 


M cr 


M B , cr 


M 


0.0 


1.573 


1.752 


1.573 


1.0 


1.494 


1.643 


1.485 


2.0 


1.390 


1.492 


1.361 



Table 3 

Gravitational (M cr ) and baryonic (Mfl lCr ) critical mass 
(see text for more details) for proto-hadronic stars at 
different entropy per baryon S/kB- M denotes the grav- 
itational mass of the cold hadronic configuration with 
the same stellar baryonic mass (Mg tCr ). Stellar masses 
are in units of the solar mass, Mq = 1.989 X 10 33 g. 

S/ks — (corresponding to a cold hadronic 
star) i 1 and 2. In the same table, we also 
report the value of the gravitational mass M. 
of the cold hadronic star with baryonic mass 
equal to Me iCr . This configuration is stable 
(r = oo) with respect to Q*-matter nucleation 
in the case S/ks = 2, and it is essentially 
stable (having a nucleation time enormously 
larger than the age of the universe) in the case 
S/ks = 1- Note that these numbers are model 
dependent and, therefore, one must take them 
just as indicative values. A careful analysis is 
beyond the scope of the present work and it 
will be addressed in a future work. 

In summary, in this work we have studied the 
quark deconfinement phase transition in hot (3- 
stable hadronic matter, and we have explored 
some of its consequences for the physics of neu- 
tron stars at birth. Our main finding is that 
proto-hadronic stars with a mass lower that the 
critical value M cr could survive the early stages 
of their evolution without decaying to a quark 
star. However, the prompt formation of a crit- 
ical size drop of quark matter could take place 
when M > M cr . These proto-hadronic stars 
evolve to cold and deleptonized quark stars, or 
collapse to a black holes. Finally, if quark mat- 
ter nucleation occurs during post-bounce stage 



2 Notice that in ref. 1 2 3 4 5 6 the critical mass for 
cold (T = 0) metastable hadronic stars has been defined 
as the value of the gravitational mass for which the 
quantum nucleation time is equal to one year: M cr (T = 
0) = M(r q = 1 yr). It is worth recalling that the nu- 
cleation time T q is an extremely steep function of the 
hadronic star mass 1 1 121 . therefore the exact value of r q 
chosen in the definition of M cr (T = 0) is not crucial 
(one must take a "reasonable small" value of t 9 , much 
shorter that the age of young pulsars as the Crab pul- 
sar). We have verified that changing r q from 1 yr to 
10 3 s modifies M cr (T = 0) by ~ 0.02%. On the other 
hand, the nucleation time r = min(T q ,T t jj) entering in 
the definition of the critical mass of proto-hadronic stars 
M cr (S) must be comparable to the proto-hadronic star 
cooling time t coo i. 
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of core-collapse supernova, then the quark de- 
confinement phase transition could trigger a de- 
layed supernova explosion characterised by a 
peculiar neutrino signal [50] . 
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